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Abstract 

Gamow- Teller strengths for selected nuclei in the iron region (A ~ 56) have 
been investigated via shell-model Monte Carlo calculations with realistic in- 
teractions in the complete fp basis. Results for all cases show significant 
quenching relative to single-particle estimates, in quantitative agreement with 
(n,p) data. The J = 1, T = residual interaction and the fi/2~f5/2 spin-orbit 
splitting are shown to play major roles in the quenching mechanism. Calcu- 
lated B(E2,2f — > Of) values are in fair agreement with experiment using 

effective charges of e p = 1.1 e and e n = 0.1 e. 

PACS numbers: 21.60.Cs, 21.60.Ka, 27.40,+z, 25.40.Kv 
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Gamow- Teller (GT) transitions from the ground states of medium and heavy nuclei are 
important in both pre-supernova neutronization |L| and double /3-decay 0. These transitions 
can be investigated experimentally through electron capture//3-decay studies and through 
intermediate-energy (p, n) and (n, p) reactions. A long-standing puzzle has been that exper- 
iments systematically yield total GT strengths several times smaller than single-particle 
estimates. 

A number of explanations have been proposed for the observed quenching of the GT 
strength, including a depletion of strength by A-excitations at about 300 MeV [|J and a 
reduction of the axial coupling constant to Qa — 1 in nuclear matter ||. A less exotic cause 
would be the presence of multi-particle, multi-hole configurations, and in light nuclei, where 
complete Ohu shell model calculations are possible, the GT strength generally decreases 
as the model space is enlarged ||. However, an exponentially increasing computational 
difficulty limits similar studies in nuclei in mid-/p-shell to the 2p-2h level, where calculations 
show a quenching only about half of that observed experimentally. 

While full Ohu calculations apparently do not recover the complete quenching of the 
GT-strength for s<i-shell nuclei 0, there are phenomenological indications that a complete 
treatment of the fp-shell is both necessary and sufficient to describe GT quenching in the 
iron region @, and a first calculation of 54 Fe in such a basis with a realistic interaction 
resulted in a quenching significantly larger than truncated estimates and comparable to 
experiment ||. The purpose of this Letter is to present calculations for several nuclei in 
the iron region ( 54 Cr, 55 Mn, 54,56 Fe, and 56,58 Ni) aimed at exploring the universality of this 
quenching and its dependence upon the effective interaction. 

Our calculations are performed in the complete set of 0/7/2,5/2^1^3/2,1/2 configurations 
using the Monte Carlo methods described in Refs. As in previous truncated shell- 

model calculations of iron region nuclei B, we have used for the bulk of our studies the 



Brown-Richter Hamiltonian |10[ fitted to lower /p-shell nuclei, but we have also calculated 
54 Fe with the Kuo-Brown interaction [IT]. Calculations for 54,56 Fe, 56 > 58 Ni, and 54 Cr were 



performed at = 2 MeV 1 and 32 time slices (both of which were checked to be sufficiently 



large to guarantee cooling to the ground state), while our 55 Mn results are at (5 — 1 MeV -1 
and 16 time slices, evaluated using a weighted average of a 54 Cr ensemble. Each calculation 
involved some 3300 Monte Carlo samples at each of six values of the coupling constant g || 
equally spaced between —1 and 0; extrapolation to the physical case (g — 1) was done by 
the method described in Ref. J8|. 

In Table | we give results for selected static observables. When the calculated bind- 
ing energies are corrected for the Coulomb energy using the semi-empirical formula E c = 
0.717(Z 2 /A 1 / 3 )(1 - 1.69/A 2 / 3 ) MeV @, all nuclei are overbound by about 1.5-2 MeV ex- 
cept 56 Ni (~ 3 MeV) and 58 Ni (~ MeV). The B(E2, 2| Of) values calculated with bare 
nucleon charges are typically within a factor of 2 of the experimental values if we assume 
that this transition saturates the total Ohoo strength. The agreement is improved if core 
polarization effects are simulated by introducing effective charges that account for coupling 
to 2 + configurations involving other major shells. We find that adopting e p = 1.1 e and 
e„ = 0.1 e for protons and neutrons, respectively, reproduces the observed B(E2) values for 
54 Fe, 56 Ni, and 58 Ni, while for the open-shell nuclei 54 Cr and 56 Fe slightly larger effective 
charges (e p ~ 1.2 e, e n ~ 0.2 e) are needed, indicating a larger core polarization in the latter 
two than in the magic and semi-magic nuclei of the A « 56 region. These effective charges 
are significantly smaller than those required in truncated shell model calculations, indicating 
that the complete Ohu model space contains significant correlations absent in smaller model 
spaces. 

In Table m we list results for the Gamow- Teller strengths _B(GT + ) = ((a-r + ) 2 ). In general, 
the GT + strength is quenched significantly relative to the single-particle value, and is in good 
agreement with experiment in all cases except 58 Ni. Since the GT operators ar + induce 
only Oku transitions, we do not expect a significant renormalization of the operators. Note 
that the quenching factors vary significantly from one nucleus to another and are not well 
approximated by a common constant value, as is conventional in astrophysical applications 
|T3|j . Compared to restricted {2p-2h) shell model approaches ||, our full /p-shell calculation 
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recovers about twice the quenching for both 54 Fe and 56 Ni. 

The energy centroids of the GT± strengths are also listed in Table [Tl] and can be com- 
pared with the experimental results from (n,p) and (p,n) reactions ||. Since we use 
isospin-invariant Hamiltonians, analogue states are exactly degenerate and the calculated 
GT_ centroid is the energy of the GT resonance relative to the isobaric analogue state, 
AE = Eqt_ — Eias- Experimental values for the latter are well parametrized by |14 



AE = (6.8 - 27.9 ■ (N - Z)/A) MeV. Although (N - Z)/A in our case is outside the range 
studied by Nakayama et al. ||14j| , we find that the MC results in all of our cases agree with 
this parametrization to within 0.4 MeV, except for the double-magic nucleus 56 Ni where 
shell-closure effects apparently not considered in the empirical parametrization might be 
important. 

To compare the GT± centroids to the data, we calculate the excitation energy of the 
resonance E x in the daughter nucleus either by using known analogue states or from the 
experimental mass difference corrected by the Coulomb energy. Although the calculated 
centroids for the (n,p) reactions (GT + ) are systematically 1.5 or 2 MeV too low, they 
generally track the position of the GT resonance in the daughter nucleus well for all nuclei. 
In particular, the calculation places the GT strength for the odd-Z target 55 Mn at a higher 
excitation energy E x than in the neighboring nuclei, in accord with experiment. A similarly 
high excitation energy is also observed for 51 V and 59 Co, which aside from 55 Mn, are the 
only odd-Z, mid-/p-shell nuclei experimentally studied thus far, suggesting an odd-even 
dependence of the excitation energy centroid. This might have an important consequence 
for the measurement of the S(GT + ) strength in (n, p) reactions, which can determine the 
strength reliably only up to daughter excitation energies of about 8 MeV. Given the odd-even 
dependence, experiments with odd-Z targets are therefore likely to "miss" a relatively larger 
fraction of the total strength with E x > 8 MeV, an effect that must be taken into account 
if one wants to compare total strengths 0. The odd-even dependence should also be quite 
important in astrophysical applications like pre-supernova calculations, but apparently has 
been neglected to date. 



It is important to determine the causes of the large GT quenching. There is strong 
evidence that neutron-proton correlations in the ground state are the source of the quenching, 



e.g., Ref. JT5[. Engel et al. have suggested that the GT quenching is particularly sensitive to 
the T = 0, J = 1 matrix element | I6[| . To investigate this point, we show in Fig. [1] B(GT + ) 
for 54 Fe in a calculation where all T = 0, J = 1 particle-particle matrix elements have been 
scaled by < g pp < 2 (with g pp = 1 being the physical value). We show the results for g = 
only, but we expect the extrapolated results at g = 1 to exhibit a similar behavior. There is 
a high sensitivity to g pp , and, upon comparing the result at g pp = (all T = 0, J = 1 matrix 
elements vanishing) to that at g pp = 1, we see that this component of the interaction causes 
about half of the total quenching, in good agreement with the phenomenological prediction 
of Ref. |16|] . In a further calculation, we switched off allT = matrix elements and recovered 
about 90% of the single-particle estimate for B(GT + ). Hino et al. have argued that the 
increase of ground state correlations is accompanied by an increase of (Jp) and (J^) as 
proton-neutron states with J / will be admixed into the ground state by pn-correlations 
[T7J . The increase of (J pn ) with increasing g pp shown in Fig. 1 supports this argument. 

The GT quenching is also expected to be very sensitive to the j '7/2-/5/2 spin-orbit split- 
ting. To investigate this effect, we have modified the BR interaction by decreasing the 
/5/2 single-particle energy arbitrarily by 2 MeV, thus reducing £5/2 — £7/2 from 6.49 MeV 
to 4.49 MeV (this change defines the "MBR" interaction). B(GT + ) values calculated with 



this force are shown in Table |TJ, where we also list effective spin-orbit splittings derived 
from the difference in the energy centroids of the neutron spectral functions, (0^07/2) and 
( a 5/2<4/2)- The effective spin-orbit splitting generally tracks the single-particle energies, and 
is comparable to the expected single-particle splitting of 7.16 MeV. It is clear from these 
results that even modest changes in the single-particle energies can have large effects on 
the quenching. Lowering the / 7 / 2 — /s/ 2 energy splitting decreases the total binding energy 
(which is overestimated by BR) significantly. However, it also decreases the centroids of the 
GT strengths which were already too low. Thus, the overbinding and the underestimation 
of the GT + excitation energy cannot be solved simultaneously by a shift of the j '7/2-/5/2 
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energy splitting. 

We have also calculated 54 Fe with the KB interaction. We find a GT strength about 
a factor of two smaller than for the BR interaction, indicating a strong, unanticipated 
sensitivity of the GT quenching to the effective interaction. A similar dependence is also 
observed in the B(E2) values, for which the KB interaction predicts a value 50% larger 
than the BR interaction and in excess of the data. As the single particle energies of the 
two interactions are nearly equal, the observed differences must stem from differences in the 
two-body interaction. 

In summary, our complete /p-basis calculations of nuclei with A ~ 56 demonstrate 
that substantial quenching of the GT + strength is a general phenomenon for all of the 
nuclei calculated and that restricted shell model calculations miss a significant part of this 
quenching. Our calculations with the BR interaction give a quenching comparable to that 
observed in (n,p) reactions; a calculation of 54 Fe with the KB interaction gives an even 
greater quenching. Our results show that np correlations are a major contributor to this 
quenching, as is a small / 7/2-/5/2 

spin-orbit splitting. The binding energies and B(E2, 2^ — ► 
0i) values calculated with the BR interaction are in reasonable agreement with experiment. 
In particular, the latter require significantly smaller effective charges than do previous, more 
restricted calculations. 

Of course, all of our results depend upon the Hamiltonian used. The BR interaction 
has been fitted to the lower fp-shell, and is probably not optimal for the present studies. 
In fact, we find some systematic shortcomings as the BR interaction tends to overbinding 
while placing the GT centroids at somewhat too low an excitation energy. Given these 
deficiencies and the strong sensitivity of the GT + strength to the effective interaction, it 
is premature to draw definite conclusions from the apparent good agreement between the 
experimental and calculated quenching of the GT + strength. A definite conclusion as to 
whether the quenching can be totally recovered in complete /p-basis calculations must wait 
until improved effective interactions are available. The new computational techniques we 
have used in this work would allow such an interaction to be determined in the complete fp 
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or even fpgg/2 basis. 
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FIGURES 

FIG. 1. The Gamow- Teller strength B(GT + ) and the expectation values of (J 2 ) for protons 
(full circles) and neutrons (open circles) as a function of the scaled strength of the T = 0, J = 1 
interaction, g pp , defined in the text. The calculation has been performed for 54 Fe at (3 = 2 MeV -1 
using the BR interaction and for the coupling constant g = 0. 
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TABLES 



TABLE I. Static observables calculated for selected nuclei with A ~ 56. 



Nucleus 


54p e 


54p e 


56 Ni 


54 Cr 


55 Mn 


56 Fe 


58 Ni 


Force 


KB a 


BR 


BR 


BR 


BR 


BR 


BR 


-(H) (MC) b 


— 


131.4 ±0.4 


145.2 ±0.6 


134.4 ±0.5 


141.8 ±1.0 


151.7±0.6 


164.4 ±0.7 


(exp) b 


129.7 


129.7 


141.9 


132.0 


140.0 


150.2 


164.4 


(Q 2 ) (fm 4 ) c 


2560 ±83 


1482 ±84 


1572 ±13 


1408 ±64 


1447 ±55 


1819 ±91 


1674 ±18 


{QD (fm 4 ) c 


368 ±31 


381 ±34 


380 ±3 


424 ±36 


384 ±24 


416 ±41 


520 ±7 


(Ql) (fm 4 ) c 


718 ±29 


478 ±5 


487 ±3 


420 ±23 




542 ±32 


528 ±6 


{QD (fm 4 ) c 


749 ±26 


454 ±30 


487 ±3 


496 ±31 




582 ±37 


569 ±7 


B(E2) (MC) d 


144 ±6 


96 ±1 


98 ±1 


84 ±5 




108 ±6 


106 ±1 




199 ±9 


129 ±1 


132 ±1 


114 ±6 




148 ±9 


142 ±2 


(exp) d 


124 ±10 


124 ±10 


120 ±10 


174 ±8 




196 ±8 


139 ±4 


(Mi 2 )Hr 


11.9 ±0.4 


14.1 ±0.4 


17.7±0.2 


13.2 ±0.4 


21.6 ±0.4 


14.7 ±0.4 


17.5 ±0.4 


a The KB interaction has 


an undefined energy shift, 


so that (H\ 


> is meaning 


dess. 





b Binding energies (in MeV) relative to 40 Ca. 

°Q p and Q n are the proton and neutron mass quadrupole moments, calculated with oscillator 
length b = 1.97(A/56) 1 / 6 fm, while Q = Q P + Q n and Q v = Q p - Q n . 

d B(E2, 2± — > 0^) in units of e 2 fm 4 calculated with the bare charges (upper row) and effective 
charges e p = Lie and e n = O.le (lower row) and assuming that this transition saturates (Q 2 ). 
e Ml is the magnetic moment operator assuming free nucleon (/-factors. 
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TABLE II. Gamow- Teller strengths for various nuclei with A ~ 56. Upper section: the Monte 
Carlo (MC) results for the total GT + strengths are compared with the single-particle values and 
the experimental strengths extracted from (n,p) data. The lower two sections give the calculated 
energy centroids of the GT strengths relative to the target ground state (Egt) an d in the daughter 
nucleus (E x ). The experimental data are from ||. The quoted values for Eqt + assume isospin 
conservation. 



Nucleus 


54 Fe 


54 Fe 


56 Ni 


54 Cr 


EC, , 

55 Mn 


56 Fe 


58 Ni 


Force 


KB 


BR 


BR 


BR 


BR 


BR 


BR 


(GT* ) (sp) 


10.3 


10.3 


13.7 


6.9 


8.6 


10.3 


13.7 


(GT^) (MC) 


2.2 ±0.3 


4.3 ±0.2 


7.4 ±0.3 


1.4 ±0.1 


2.2±0.2 


2.73 ±0.04 


5.6 ±0.3 


(GT= ) (exp) 


3.5±0.7 a 


3.5±0.7 a 






1.72±0.2 b 


2.85 ±0.3 b 


3.76 ±0.4 b 


E G t + (MC) 


9.2 ±0.2 


9.7 ±0.2 


8.9 ±0.3 




12.4 ±1.6 


11.0 ±0.2 


9.8 ±0.2 


E x (MC) 


0.7±0.2 


1.2 ±0.2 


2.7 




2.4±1.6 


-0.4 ±0.2 


1.2 ±0.2 


E x (exp) 


2.8 


2.8 






4.1 ±0.5 


2.7 ±0.5 


3.5 ±0.5 


AE (MC) 


5.2 ±0.2 


6.1 ±0.2 


8.9 ±0.3 






4.4 ±0.1 


6.1 ±0.2 


(sys) c 


5.8 


5.8 


6.8 


3.7 


4.3 


4.8 


5.9 


E x (MC) 


5.2 ±0.2 


6.1 ±0.2 








7.8 ±0.1 


6.15 ±0.2 


E x (exp) 


8.2 


8.2 













a Summed to an excitation energy of 9 MeV. 
b Summed to an excitation energy of 8 MeV. 



"The systematics of Ref. |14j for ^gt_ — -Etas 
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TABLE III. Selected properties of 54 Fe and 58 Ni calculated at (5 = 2 with the Brown-Richter 
(BR) interaction of Ref. [|l(| and its modification (MBR), where the _f 5 / 2 single-particle energy has 
been lowered by 2 MeV. 



Nucleus 


^e 


54-cu 
re 


58 at; 
JNi 


JNl 


Force 




Mrjrt 


dK 


IvioK 






An +n 

Li.^iX) xu.uo 


K ft _i_n Q 


9 qc _i_n f)A 




9.7±0.2 


8.2±0.2 


9.8 ±0.2 


7.2 ±1.0 


-(H) 


131.4 ±0.4 


126.4 ±0.5 


164.4 ±0.7 


159.8 ±0.5 




11.77 ±0.08 


12.36 ±0.03 


11.73 ±0.09 


10.40 ±0.03 




-3.45 ±0.08 


-5.54 ±0.08 


-3.36 ±0.16 


-4.18 ±0.09 




8.31 ±0.11 


6.82 ±0.08 


8.37 ±0.18 


6.22 ±0.09 


a Energy centroid 


in MeV of the 


neutron / 7 / 2 strength, 


obtained from the 


response function 



(4 /2 (r)a 7/2 ). 

b Energy centroid of the neutron / 5 / 2 strength, obtained from (a 5 / 2 (r)agy 2 ). 
c Effective spin-orbit splitting, E 7 / 2 + E 5 / 2 - 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9402006vl 



